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In this paper we introduce a critical curve separating the asymp- 
totic behavior of the moments of the symbiotic branching model, 
introduced by Etheridge and Fleischmann [Stochastic Process. Appl. 
114 (2004) 127-160] into two regimes. Using arguments based on 
two different dualities and a classical result of Spitzer [Trans. Amer. 
Math. Soc. 87 (1958) 187-197] on the exit-time of a planar Brownian 
motion from a wedge, we prove that the parameter governing the 
model provides regimes of bounded and exponentially growing mo- 
ments separated by subexponential growth. The moments turn out to 
be closely linked to the limiting distribution as time tends to infinity. 
The limiting distribution can be derived by a self-duality argument 
extending a result of Dawson and Perkins [Ann. Probab. 26 (1998) 
1088-1138] for the mutually catalytic branching model. 

As an application, we show how a bound on the 35th moment 
improves the result of Etheridge and Fleischmann [Stochastic Process. 
Appl. 114 (2004) 127-160] on the speed of the propagation of the 
interface of the symbiotic branching model. 



1. Introduction. In 2004, Etheridge and Fleischmann [8] introduced a 
stochastic spatial model of two interacting populations known as the sym- 
biotic branching model, parametrized by a parameter g& [—1)1] governing 
the correlation between the two driving noises. The model can be considered 
in three different spatial setups which we now explain. 
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First, the continuous-space symbiotic branching model is given by the 
system of stochastic partial differential equations 

d 



(1.1) cSbm(£»,k; 



Uq,Vq 



: < 



1 



dt 
d 



ut{x) = -Aut{x) + -s/ KUt{x)vt{x) dWl{x) 



Q^Vt{x) = -Avt{x) + yjKUt{x)vt{x)dW'l{x), 

uq{x) > 0, 
. vo{x) > 0, 



X E . 
X £ 



where A denotes the Laplace operator and k > is a fixed constant known as 
the branching rate. W = (W^, W'^) is a pair of correlated standard Gaussian 
white noises on R_|_ x M with correlation g G [—1,1], that is, the unique 
Gaussian process with covariance structure 



(1.2) 
(1.3) 
(1.4) 



nwl{A,)Wl{A2)] 
E[Wl{A,)Wl{A2)] 
nwl{A^)Wl{A2)] 



{tiAt2)i{AinA2), 
{tiAt2)i{AinA2), 
g{h At2)i{Air\A2) 



where i denotes Lebesgue measure, ^i, G B{M.) and ti,t2 > 0. Note that 
we work with a white noise W in the sense of Walsh [25] . Solutions of this 
model have been considered rigorously in the framework of the corresponding 
martingale problem in Theorem 4 of [8], which states that, under suitable 
conditions on the initial conditions uq{-),vq{-), a solution exists for all gG 
[—1,1]. The martingale problem is well posed for all ^ G [—1,1), which implies 
the strong Markov property except in the boundary case g = l. 

For a discrete spatial version we consider the system of interacting diffu- 
sions on Z'^, with values in M>o, defined by the coupled stochastic differential 
equations 



(1.5) dSBM(^,K; 



uo ,vo 



: < 



where now {i?^(z), i?^(i)}jg 
covariances given by 



dut{i) = Aut{i) dt + KUt{i)vt{i) dBl{i), 

dvt{i) = Avt{i) dt + nut{i)vt{i) dB'^{i), 
no(i)>0, ieZ"^, 
uo(i)>0, ieZ'^, 

d is a family of standard Brownian motions with 



(1.6) 



[B-{t),B^{j)]t 



i= j and n^m, 
i = j and n = m, 
otherwise. 



In the discrete case, A denotes the discrete Laplacian 

1 



Autii)= ^(utik) - ut{i)) 



|fc-i|=i 
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Note that in this paper we denote by [N.,M.]t the cross- variation of two 
martingales iV, M . This is to avoid confusion with (/, g) which will be defined 
to be the sum (resp., integral) of the product of / and g. 

Finally, the nonspatial symbiotic branching model is defined by the stochas- 
tic differential equations 



Again, the noises are correlated with [B^^B'^\t = Qi- This simple toy-model 
(see also [19] and [6]) can be analyzed quite simply and will be used to prove 
properties of the spatial models. 

Convention 1.1. Prom time to time we skip the dependence on q^k, uq 
and vq if there is no ambiguity. Solutions of cSBM, SBM and dSBM for 
d < 2 are called symbiotic branching processes in the recurrent case whereas 
solutions of dSBM for d>2> are called symbiotic branching processes in the 
transient case. 

Interestingly, symbiotic branching models include well-known spatial mod- 
els from different branches of probability theory. In the discrete spatial case 
(and analogously in continuous-space) interacting diffusions of the type 



have been studied extensively in the literature. Some important examples 
are the following: 

Example 1. The stepping stone model from mathematical genetics: 



Example 2. The parabolic Anderson model (with Brownian potential) 
from mathematical physics: f{x) = . 

Example 3. The super random walk from probability theory: f{x) = x. 

For the super random walk, k is the branching rate which in this case is 
time-space independent. In [7], a two-type model based on two super random 
walks with time-space dependent branching was introduced. The branching 
rate for one species is proportional to the value of the other species. More 
precisely, the authors considered 




(1.7) 



dwt{i) = Awt{i) dt + -s/ Kf{wt{i)) dBt{i) 



f{x) = x{l-x). 



dut{i) = Aut{i) dt + sj KUt{i)vt{i) dB}{i) 



dvt{i) = Avt{i) dt + ^ K.ut{i)vt{i) dBf{i), 
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where now {B^{i),B^{i)}j^i^^d is a family of independent standard Brownian 
motions. Solutions are called mutually catalytic branching processes. In the 
following years, properties of this model were well studied (see, e.g., [3] and 
[2]). The corresponding continuous-space version was also treated in [7]. 

For correlation Q = 0, solutions of the symbiotic branching model are 
obviously solutions of the mutually catalytic branching model. The case g = 
— 1 with the additional assumption uo + vo = l corresponds to the stepping 
stone model. To see this, observe that in the perfectly negatively correlated 
case B^{i) = —B'^{i) which implies that the sum u + v solves a discrete heat 
equation and with the further assumption uq + vq = 1 stays constant for all 
time. Hence, for all t > 0, u{t, •) = 1 — v (t, •), which shows that n is a solution 
of the stepping stone model with initial condition uq and v is a solution 
with initial condition vq. Finally, suppose w \s a, solution of the parabolic 
Anderson model, then, for ^=1, the pair {u,v) := {w,w) is a solution of the 
symbiotic branching model with initial conditions uq = vq = wq. 

The purpose of this and the accompanying paper [1] is to understand the 
nature of the symbiotic branching model better. How does the model depend 
on the correlation g? Are properties of the extremal cases g € { — 1,0, 1} in- 
herited by some parts of the parameter space? Since the longtime behavior 
of the super random walk, stepping stone model, mutually catalytic branch- 
ing model and parabolic Anderson model is very different, one might guess 
that the parameter space [—1, 1] can be divided into disjoint subsets corre- 
sponding to different regimes. 

The focus of [1] is second moment properties. In the discrete setting, but 
with a more general setup, growth of second moments is analyzed in detail. 
A moment duality is used to reduce the problem to moment generating 
functions and Laplace transforms of local times of discrete-space Markov 
processes. A precise analysis of those is used to derive intermittency and 
aging results which show that different regimes occur for g <0, g = and 
g>0. 

In contrast to [1], the present paper is not restricted to second moment 
properties. The aim is to understand the pathwise behavior of symbiotic 
branching processes better. 

Remark 1.2. In this paper, we restrict ourselves to the simplest setups 
which already provide the full variety of results. For the discrete spatial 
model we thus restrict ourselves to the discrete Laplacian instead of al- 
lowing more general transitions. This is not necessary; see [7] or [2] for a 
construction of solutions and main properties for more general underlying 
migration mechanisms in the case = 0. Furthermore, we mainly restrict 
ourselves to homogeneous initial conditions and remark where results hold 
more generally. Here, for nonnegative real numbers we denote by u the con- 
stant functions u{-)=u. 
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The paper is organized as follows: our main results are presented in Section 
2. Before proving the results, we collect basic properties of the symbiotic 
branching models and discuss the dualities that we need. This is carried out 
in Section 3. The final sections are devoted to the proofs. In Section 4, proofs 
of the longtime convergence in law are given, and in Section 5 we discuss 
the longtime behavior of moments. Finally, in Section 6 we show how to use 
the results of Section 5 to strengthen the main result of [8]. 

2. Results. Before stating the main results, we briefly recall from [8] 
that the state space of cSBM is given by pairs of tempered functions, that 
is, pairs of functions contained in 

-^tem = : IK — ^ IK>Oi li™ exists and ||ii(7!)A||oo < oo VA < [•, 

where (j)x{x) = e'^'^' , and we think of Mtcm as being topologized by the metric 
given in [8], equation (13), yielding a Polish space. 

The state space for dSBM is similar. It was not discussed in [8] and so we 
present details in Section 3. 

2.1. Convergence in law. We begin with a result, generalizing Theorem 
1.5 of [7], on the longtime behavior of the laws of symbiotic branching pro- 
cesses in the recurrent case. 

Proposition 2.1. Suppose {ut,vt) is a spatial symbiotic branching pro- 
cess in the recurrent case with g G ( — 1,1), n > and initial conditions 
uq = u,vo = v. Let and be two Brownian motions with covariance 

[B^,B% = et, t>0, 

and initial conditions Bq = u,Bq = v. Further, let 

T = mf{t > : BI Bf = 0} 

be the first exit time of the correlated Brownian motions B^,B^ from the 
upper right quadrant. Then, weakly in M^^^, 

F-^^[{ut,Vt)G-]^P^'^m,B',)€-] 

as t ^ oo. Here, {BI^B"^) denotes the pair of constant functions on M, re- 
spectively, Z'^ (d=l,2) taking the values of the stopped Brownian motions 
{Bl,Bl). 

In particular, the proposition shows ultimate extinction of one species in 
law. 
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Remark 2.2. For simplicity, Proposition 2.1 is formulated for constant 
initial conditions even though the result holds more generally. Theorem 1.5 of 
[7] (the case g = 0) was extended in [4] to nondeterministic initial conditions: 
for fixed u,v'>0 let M-u,v be the set of probability measures v on M^^^ such 
that 

(2.1) sup / (a^(x) + b^{x)) du{a, b) < oo 

x&rJ 

and 

(2.2) lim / [{Pta{x) - uf + (Pj6(x) - vf] du{a, b) = for all x G M. 

t-^oo J 

Here, (Pt) denotes the transition semigroup of Brownian motion (the defi- 
nition for the discrete case is similar). The proof of [4] can also be applied 
to g^O and, thus, Proposition 2.1 holds in the same way for initial distri- 
butions v S J^u,v 

The restriction to g £ (—1, 1) arises from our method of proof which ex- 
ploits a self-duality of the process which gives no information for g€ {—1,1}. 
Let us briefly discuss the behavior of the limiting distributions in the bound- 
ary cases € {—1,1} which are well known in the literature and fit neatly 
into our result. First, suppose (wt) is a solution of the stepping stone model 
(see Example 1) and wo = w£ [0, 1]. It was proved in [20] that 

(2.3) C'^iwt)''^w6i + {l-w)So, 

where 5± (resp., 6q) denotes the Dirac distribution concentrated on the con- 
stant function 1 (resp., 0). This can be reformulated in terms of perfectly 
anti-correlated Brownian motions {B^,B'^) as before: for q= —1, the pair 
{B^,B'^) takes values only on the straight line connecting (0,1) and (1,0), 
and stops at the boundaries. Hence, the law of {Bl,B^) is a mixture of ^(o,i) 
and (5(1^0) the probability of hitting (1,0) is equal to the probability of 
a one-dimensional Brownian motion started in w £ [0, 1] hitting 1 before 0, 
which is w, and hence matches (2.3). Second, let (wt) be a solution of the 
parabolic Anderson model with Brownian potential (see Example 2) and 
constant initial condition wq = w >0. In [21] it was shown that 

(2.4) C'^{wt)'^6o. 

As discussed above, when viewed as a symbiotic branching process with 
g= I, this implies 

(2.5) £-'-(nt,^;t)*^°°5o,o. 

From the viewpoint of two perfectly positive-correlated Brownian motions, 
we obtain the same result since they simple move on the diagonal dissecting 
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the upper right quadrant until they eventually get absorbed in the origin, 
that is, {BI.BI) = (0,0) almost surely. 

To summarize, we have seen that the weak longtime behavior (in the re- 
current case) of the classical models connected to symbiotic branching is 
appropriately described by correlated Brownian motions hitting the bound- 
ary of the upper right quadrant. 

2.2. Nonalmost-sure behavior. In contrast to extinction in law, the almost- 
sure behavior is very different. In the recurrent case for the mutually cat- 
alytic branching model. Cox and Klenke [3] showed that, almost surely, there 
is no longtime local extinction of any type, but in fact the locally predom- 
inant type changes infinitely often. It is not hard to see that the same is 
true for symbiotic branching with q ^ (—1,1). We do not give a proof since 
it follows from Proposition 2.1 along the same lines as in [3]. 

Proposition 2.3. Let g £ (—1,1), k > and suppose {ut,vt) is a spatial 
symbiotic branching process in the recurrent case with initial distribution 
uq = u,tio = V. Then, for all {u',v') € {{x,0):x € M>o} U {{0,y):y G M>o} 
and K cM. bounded, 



liminf sup||(uj(x),'yt(x)) — (u',f')|| = 

L t->oo 



1, 



respectively, for K dU^ bounded. 



liminf sup||(tij(/c),tij(A;)) — {u' ,v' 
L t^oo fcg^ 



Again, as in Remark 2.2, the result holds for random initial conditions of 
the class Mu,v Note that Proposition 2.3 depends strongly on the spatial 
structure since in the nonspatial model almost sure convergence holds (see 
Proposition 4.4). 

2.3. Longtime behavior of moments. In [1] the second moments of symbi- 
otic branching processes are analyzed. This particular case admits a detailed 
study since a moment duality (see Lemma 3.3) has a particularly simple 
structure which allows one to reduce the study of the moments to that of 
moment generating functions and Laplace transforms of local times. Here 
we are interested in the behavior of moments as t tends to infinity. The 
two available dualities (self-duality and moment duality) are combined in 
two steps. First, a self-duality argument combined with an equivalence be- 
tween bounded moments of the exit time distribution and of the exit point 
distribution for correlated Brownian motions stopped on exiting the first 
quadrant is used to understand the effect of q. It turns out that for any 
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Fig. 1. The critical curve p{q), g £ ( — 1,1). 



p > 1 there are critical values, independent of k, dividing regimes in which 
the moments E^'^[uf], ¥}^'^[ut{ky] and ¥?-^^[ut{xy] are bounded in t or grow 
to infinity. Second, for p G N, a perturbation argument combined with the 
first step and a moment duality is used to analyze the growth to infinity in 
more detail. 

The following critical curve captures the effect of q. Note that the defi- 
nition is independent of k which will become important in the second step. 



Definition 2.4. We define the critical curve of symbiotic branching 
models to be the real- valued function p: (— 1, 1) — t- M"*", given by 



7r/2 + arctan(£»/(v^l - q^)) ' 
Its inverse will be denoted by q{p) for p > 1. 

The critical curve is plotted in Figure 1. Here, f?(35) and ^(2) are marked. 
Thirty-fifth moments are the key for the improved wavespeed result below 
and the special case q{2) = is discussed in [1]. We will see in Section 5 
that this curve is closely connected with the exit distribution of {Bl,B^) 
from the upper right quadrant which appeared in Proposition 2.1 above. 
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The first main theorem states that the critical curve separates two regimes 
(independently of k): that of bounded moments and that of unbounded 
moments. 

Theorem 2.5. Suppose {ut,vt) is a symbiotic branching process (either 
nonspatial, continuous space or discrete space in arbitrary dimension) with 
initial conditions uq=vq = '\-. If g £ {—1,1), then, for any k> 0, the follow- 
ing hold for p> 1: 

(i) In the recurrent case, 

q<q{p) ^ E^'^[u^],E^^^[utik)P] and E^^^[ut{x)P] are bounded int. 

(ii) In the transient case, 



Due to symmetry the same holds for K^'^lvf], K^''^[vt{k)'^] and ¥,^''^[vt{xy]. 

Note that the theorem provides information about all positive real mo- 
ments, not just integer moments. In the area below the critical curve in 
Figure 1, the moments remain bounded. On and above the critical curve, in 
the recurrent case, the moments grow to infinity. 

Remark 2.6. For q = —1 the curve could be extended with p{—l) = oo. 
In terms of the previous theorem this makes sense since for q = —1, symbiotic 
branching processes with initial conditions uq = vq = 1 are bounded by 2. 
This is justified by a simple observation: for initial conditions uq = vo = 1/2 
symbiotic branching processes with g = —1 are solutions of the stepping 
stone model and, hence, bounded by 1. Uniqueness in law of solutions implies 
that solutions [ut,vt) with initial conditions {cuq,cvq) are equal in law to 
solutions c times solutions with initial conditions {uq,vq). 

With this first understanding of the effect of g on moments, we may dis- 
cuss integer moments for the discrete-space model in more detail. Let us first 
recall some known results for solutions {wt) of the parabolic Anderson model 
(see Example 2) where only the parameter k appears. Using Ito's lemma, one 
sees that m{t,ki, . . . ,kn) :=K^[wt{ki) ■ ■ -wtikn)] solves the (discrete-space) 
partial differential equation 

d 

—m{t,ki, ...,kn) = Am{t,ki, ...,kn) + V{ki, . . . , fc„)m(t, /ci, . . .,kn) 
with homogeneous initial conditions. Here, the potential V is given by 



Q< q{p) =^ ^'^'^[ut{kY] is bounded in t. 




l<i<j<n 
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Since H = —A — y is an n-particle Schrodinger operator, many properties 
are known from the physics hterature. In particular, it is well known that 
in the recurrent case (the potential is nonnegative) exponential growth of 
solutions holds for any k > 0. By contrast, in the transient case the discrete 
Laplacian requires a stronger perturbation before we see exponential growth. 
Intuitively from the particle picture this should be true since the potential 
V only increases solutions if particles meet, which occurs less frequently in 
the transient case. For the transient case (see, e.g., [5] or [9] for more precise 
results), there is a decreasing sequence n{n) such that 

E^[wj(/c)"] is bounded in t <^ k < Hi{n) 

and for the Lyapunov exponents 

7„(k) := lim - logE^ [wtik)"-] > k> K{n). 

t— >oo t 

These results can be proved with the n-particle path-integral representation 
in which solutions are expressed as 

where (X/ ),..., (X") are independent simple random walks started in 
ki, . . . ,kn- 

Coming back to the symbiotic branching model, we ask whether or not 
the nth Lyapunov exponents 

7n(£',«) := lim ^logE^ '\ut{kr] 

exist and in which cases 7n(f?)^) is strictly positive. As for the parabolic 
Anderson model, there is a system of partial differential equations describing 
the moments (see Proposition 16 of [8] for the continuous-space model) and 
an n-particle path-integral representation of the moments. In addition to 
the independent motion, the particles now carry a color which randomly 
changes if particles of the same color stay at the same site (see Lemma 3.3). 
With denoting collision times of particles of same color and Lf denoting 
collision times of particles of different colors, the path- integral representation 
of moments reads 

Ei'i[nt(A;)"]=E[e'^(-^^+^-^^)]. 

This representation is more involved than the path-integral representation 
for the parabolic Anderson model since, in addition to the motion of par- 
ticles, a second stochastic mechanism is included. Nonetheless, we use it to 
prove the following theorem which reveals that even in the recurrent case a 
nontrivial transition occurs. 
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Theorem 2.7. For solutions 0/ dSBM(^, fi;)i^i, in any dimension, the 
following hold for n G N, n > I.- 

(i) ln{Q-,K) exists for any Q ^[—1,1], k>0, 

(ii) 7n = for any k>0, 

(iii) for any g > g{n) there is a critical K{n) such that jniQ,/^) > ^ if 
K > K.{n). 

Combined with Theorem 2.5, parts (ii) and (iii) emphasize the "critical- 
ity" of the critical curve: for g < g{n), moments stay bounded, for g = g{n) 
moments grow subexponentially fast to infinity, and for g > g{n) moments 
grow exponentiahy fast if k is large enough. 

Remark 2.8. As discussed above, for the parabolic Anderson model 
it is natural that in the transient case perturbing the critical case does not 
immediately yield exponential growth, whereas perturbing the recurrent case 
does immediately lead to exponential growth. It is clear that in the transient 
case the gap in (iii) of Theorem 2.7 is really necessary: for small k moments 
of the parabolic Anderson model are bounded. Since moments of symbiotic 
branching are dominated by moments of the parabolic Anderson model (see 
Lemma 3.3), for small k moments are bounded for all g. 

In the case p ^ N there seems to be no reason why exponential growth 
should fail. Unfortunately, in this case there is no moment duality and hence 
the most useful tool to analyze exponential growth is not available. 

Conjecture 2.9. In the recurrent case the moment diagram for sym- 
biotic branching (Figure 1) describes the moments as follows: pairs {g,p) 
below the critical curve correspond precisely to bounded moments, pairs at 
the critical curve correspond to moments which grow subexponentially fast 
to infinity and pairs above to the critical curve correspond to exponentially 
growing moments. 

A deeper understanding of the Lyapunov exponents as functions of g, k 
remains mainly open (for an upper bound see Proposition 5.3). For second 
moments [g{2) =0] this is carried out in [1]. It is shown that exponential 
growth holds for g > and arbitrary k > in the recurrent case, whereas 
only for k, > 2/{gGooiO, 0)) in the transient case. Here Goo denotes the Green 
function of the simple random walk. The exponential (and subexponential) 
growth rates were analyzed in detail by Tauberian theorems. 

A direct application of Theorem 2.7 is so-called intermittency of solutions. 
One says a spatial system with Lyapunov exponents 7p is p-intermittent if 

> ^ 7p+i 
p p + l' 
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Intermittent systems concentrate on few peaks with extremely high intensity 
(see [10]). The results above show that as g tends to —1, solutions (at least 
for large k) are p-intermittent for p tending to infinity. This holds since 
for fixed the pth moments are bounded if {g,p) lies below the critical 
curve. Increasing p (and k if necessary) there is a first p such that the pth 
Lyapunov exponent is positive. Intermittency for higher exponents suggests 
that the effect gets weaker. This is to be expected since for = — 1 solutions 
with homogeneous initial conditions are bounded and, hence, solutions do 
not produce high peaks at all. Making this effect more precise, in particular 
combined with the effect of Proposition 2.1, is an interesting task for the 
future. 

2.4. Speed of propagation of the interface. Let us conclude with a di- 
rect application of the moment bounds. Here, we will be concerned with an 
improved upper bound on the speed of the propagation of the interface of 
continuous-space symbiotic branching processes which served to some extent 
as the motivation for this work. To explain this, we need to introduce the 
notion of the interface of continuous-space symbiotic branching processes 
introduced in [8]. 

Definition 2.10. The interface at time t of a solution {ut,vt) of the 
symbiotic branching model cSBM(^3, k)^;,^^,;, with g£ [—1,1] is defined as 

Ifct = cl{2; : ut{x)vt{x) > 0}, 

where cl{A) denotes the closure of the set A in JR. 

In particular, we will be interested in complementary Heaviside initial 
conditions 

uq{x) = IjR- (x) and vo{x) = 1^+ (x), xe R. 

The main question addressed in [8] is whether for the above initial conditions 
the so-called compact interface property holds, that is, whether the interface 
is compact at each time almost surely. This is answered affirmatively in 
Theorem 6 in [8], together with the assertion that the interface propagates 
with at most linear speed, that is, for each g G [—1, 1] there exists a constant 
c > and a finite random-time Tq so that almost surely for all T >Tq 

y IfctC[-cr, cT]. 

t<T 

Heuristically, due to the scaling property of the symbiotic branching model 
(Lemma 8 of [8]) one expects that the interface should move with a square- 
root speed. Indeed, with the help of Theorem 2.5 one can strengthen their 
result, at least for sufficiently small g, to obtain almost square-root speed. 
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Theorem 2.11. Suppose {ut,vt) is a solution of cSBM{q, with 
Q < ^(35) and k > 0. Then there is a constant C > and a finite random- 
time Tq such that almost surely 

U ifct c [-cVTiog{T), cVriog(r)] 

t<T 

for all T > Tq. 

The restriction to g < g{35) is probably not necessary and only caused by 
the technique of the proof. Though ^*(35) « —0.9958 is rather close to — 1, 
the result is interesting. It shows that sub-linear speed of propagation is not 
restricted to situations in which solutions are uniformly bounded as they 
are for q = —1. The proof is based on the proof of [8] for linear speed which 
carries over the proof of [24] for the stepping stone model to nonbounded 
processes. We are able to strengthen the result by using a better moment 
bound which is needed to circumvent the lack of uniform boundedness. 

Remark 2.12. We believe that, at least for g <0, the speed of propa- 
gation should be at most C'\/t, for some suitable constant C, that is, for 
all T greater than some T' > 0, 

\JUct^[-C'VT,C'VT]. 

t<T 

However, it seems unclear how to obtain such a refinement of Theroem 
2.11 based on our moment results and the method of [24] (resp., [8]). As 
subexponential bounds of higher moments cannot be avoided (see the proof 
of the fluctuation term estimate Lemma 6.2), our results on the behavior 
of higher moments show that at present, in light of Conjecture 2.9, one can 
only hope for stronger results for very small g. 

To overcome this limitation, new methods need to be employed. The au- 
thors think that a possible approach could be based on the scaling property 
(Lemma 8 of [8]) and recent results by Klenke and Oeler [13]. Recall that 
the scaling property states that if {ut,vt) is a solution to cSBM(£>, k)^;,^^,^, 
then 

{ut{xf,vt{xf) := iuKt{VKx),VKt{VKx)), xeRJOO, 

is a solution to cSBM(£>, K ■ k)^k (with suitably transformed initial states 

Uq ,Vq). In other words, a diffusive time-space rescaling leads to the original 
model with a suitably increased branching rate n. Klenke and Oeler [13] 
show that, at least for the mutually catalytic model in discrete space, a 
nontrivial limiting process as k — )• oo exists. This limit is called "infinite 
rate mutually catalytic branching process" (see also [11, 12] for a further 
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discusion). In particular, in Corollary 1.2 of [13] they claim that, under 
suitable assumptions, a nontrivial interface for the limiting process exists, 
which would in turn predict a square-root speed of propagation in our case. 
However, to make this approach rigorous is beyond the scope of the present 
paper. 

Remark 2.13 (Shape of the interface). Note that our results give only 
limited information about the shape of the interface. For the case q = —1, 
that is, with locally constant total population size, it is shown in [16] that 
there exists a unique stationary interface law, which may therefore be inter- 
preted as a "stationary wave" whose position fluctuates at the boundaries, 
according to [24], like a Brownian motion, hence explaining the square-root 
speed (note that for both results, suitable bounds on fourth mixed moments 
are required). However, for ^ > — 1, the population sizes of the interface are 
expected to fluctuate significantly and it seems unclear how this affects the 
shape and speed of the interface, in particular the formation of a "station- 
ary wave." The significance of fourth mixed moments might even lead to a 
phase-transition in g. This gives rise to many interesting open questions. 

3. Basic properties and duality. In this section we review the setting 
and properties of the discrete-space model, whereas for continuous-space we 
refer to [8] . Note that instead of using the state space of tempered functions 
alternatively we may use a suitable Liggett-Spitzer space. As the results are 
only presented for the discrete Laplacian this does not play a crucial role. 
For a discussion of the mutually catalytic branching model in the Liggett- 
Spitzer space see [2]. 

3.1. Basic properties. For functions /, abbreviate {f,g) = 

Efc/(^)5(^)- With cf>x{k) = e^\''\ the space of pairs of tempered sequences is 
defined by 

= {(.u, v)\u, v.Z'^^ M>o, {u, 0a), {v, (t>x) < oo VA < 0}. 
The space of continuous paths is denoted by 

f^tcm = C'(IR>05 M^cm)- 

Similarly, the space of pairs of rapidly decreasing sequences is defined by 

Mlp = {{u,v)\u,v:Z'^^R>o,{u,^x),{v,(l)x) <ocyX>0} 
and the corresponding path space by 

f]rap = C(M>o,M,ip). 

Weak solutions are defined as in [7] for g = 0. In much the same way as 
for Theorems 1.1 and 2.2 of [7], we obtain existence and the Green-function 
representation. 
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Proposition 3.1. Suppose {uq,vq) ^ M'^^^ (resp., M'^^^), £ig[-1,1] 
and K > 0. Then there is a weak solution o/dSBM(f3, k)uq,vo such that {ut,vt) G 
Qtem (resp., $7rapj and for all {(f),!!:) £ (resp., M^^^) 



(3.1) (uj» = (no: 



(3.2) (t;j,V') = (7;o,Pt^)+ V / PtsmV^M^M) dB^^{j), 

where Ptf{k) =J2j£ZdPtij,k)f{j) is the semigroup associated to the simple 
random walk. In particular, we have 

(3.3) Ut{k) = PtUoik) + V] [ pt-s{j,k)y/KUs{j)Vs{j)dBl{j), 

(3.4) vt{k) = Ptvoik) + V fpt-sU, k)^/KUs{j)vs{j) dBl{j). 

The covariation structure of the Brownian motions is given by (1.6). 

In fact, (3.1), (3.2) can be seen as the discrete-space versions of the martin- 
gale problem of Definition 3 in [8]. Further, (3.3), (3.4) are the discrete-space 
versions of the convolution form given in Corollary 20 of [8] . 

For the proofs of the longtime behavior of laws and moments, the key 
step is to transfer to the total mass processes (uj, 1), {vt, 1). To this end, in 
a similar way to Proposition 3.1, we define 

Ml = {{u, v)\u, v.'L'^^ M>o, (n, 1), (w, 1) < oo} 

and 

nF = C(M>o,M|). 

For summable initial conditions we obtain the following crucial martingale 
characterization . 

Proposition 3.2. //(no, uo) G M^,, then each solution of dSYiM{Q, k)uq^vo 
has the following properties: {ut,vt) G and {ut, 1), {vt, 1) are nonnegative, 
continuous, square-integrable martingales with square-functions 

[{u.,l)]t = [{v.,l)]t = K {Us,Vs)ds 



and 

[{U.,l),{v.,l)]t = 0K {Us,Vs)ds. 

Jo 
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We omit the proofs since they are basically standard. The only step where 
one needs to be careful is the existence proof. As usual for such models 
one first restricts the space to bounded subsets (boxes) of Z*^, where stan- 
dard Markov process theory applies. Enlarging the boxes one obtains a se- 
quence of processes which are shown to converge to a limiting process solving 
dSBM(^ , i^)uo,vo- To prove tightness of the approximating sequence, the mo- 
ments need to be bounded uniformly in the size of the boxes. Here, more 
care than for ^> = in [7] is needed. The uniform moment bound can, for 
instance, be achieved using a colored particle moment duality for each box 
similar to the one of Lemma 3.3. 

3.2. Dualities. The symbiotic branching model exhibits an exception- 
ally rich duality structure, providing powerful tools for the analysis of the 
longtime properties. 

3.2.1. Colored particle moment dual. We now recall the two-colors par- 
ticle moment-duality introduced in Section 3.1 of [8]. Since the dual Markov 
process is presented rigorously in [8] we only sketch the pathwise behavior. 
To find a suitable description of the mixed moment 

E«o,«o [ut{ki) ■ ■ ■ ut{kn)vt{kn+i) ■■■vt (^n+m)] 5 

n + m particles are located in U^. Each particle moves as a continuous-time 
simple random walk independent of all other particles. At time 0, n particles 
of color 1 are located at positions fci , . . . , fe„ and m particles of color 2 are 
located at positions /c„+i, . . . , For each pair of particles, one of the 
pair changes color when the time the two particles have spent in the same 
site, while both have same color, first exceeds an (independent) exponential 
time with parameter k. Let 

total collision time of all pairs of same colors up to time i, 
total collision time of all pairs of different colors up to time t, 
number of particles of color 1 at site a at time t, 
number of particles of color 2 at site a at time t, 

Note that since there are only n + m particles, the infinite product is actually 
a finite product and hence well defined. The following lemma is taken from 
Section 3 of [8]. 



Lt = 

I'M = 
(no,t;o)'* = 
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Lemma 3.3. Let {ut,vt) he a solution of dSJiM.{Q^K)uo,vQ, t > and 
^ G [-1, 1]. Then, for any h €Z'^,t>0, 

E""'"" [utih] ■ ■ ■ ut{K)vt{K+i) ■ ■ ■ vt{kn+m)] = E[(no, voy^e^^'^^+^'^t^ 
where the dual process behaves as explained above. 

Note that for homogeneous initial conditions uq = vq = 1, the first factor 
in the expectation of the right-hand side equals 1 . In the special case Q= 1, 
Uq = Vq = 1 Lemma 3.3 was already stated in [5], reproved in [9] and used 
to analyze the Lyapunov exponents of the parabolic Anderson model. 

For Q^l, the difficulty of the dual process is based on the two stochastic 
effects: on the one hand, one has to deal with collision times of random walks 
which were analyzed in [9]; additionally, particles have colors either 1 or 2 
which change dynamically. 

Remark 3.4. Similar dualities hold for cSBM and SBM. For continuous- 
space, the random walks are replaced by Brownian motions and the collision 
times of the random walks by collision local times of the Brownian mo- 
tions (see Section 4.1 in [8]). The simplest case is the nonspatial symbiotic 
branching model where the particles stay at the same site and local times 
are replaced by real times (see Theorem 3.2 of [19] or Proposition A5 of [6]). 

3.2.2. Self-duality. Mytnik [18] introduced a self-duality for the continuous- 
space mutually catalytic branching model to obtain uniqueness of solutions 
of the corresponding martingale problem. This can be extended to symbi- 
otic branching models for g G (—1, 1) as shown in Proposition 5 of [8]. The 
discrete-space self-duality for g = was proved in Theorem 2.4 of [7]. We 
first need more spaces of sequences: 

E = {{x,y): {x, \y\) € M^^, \y{k)\ < x{k) V/c G Z"] 

and 

E = {{x,y) £ E:x £ Mrap} D {{x,y) £ E:x has bounded support} = Ef. 

In the sequel, the space E and its subspaces will be used for (x,y) = {ut + 
vt,ut — vt). The duality function for g G (—1, 1) maps E x E to C via 

(3.5) H{u,v,u,v) =exp(— Y^l — g{u,u) +i\/l + g{v,v)). 

With this definition the generalized Mytnik duality states: 

Lemma 3.5. For g G (-1, 1), k>0, (uo,vo) g M^^^ and {uo,vo) G M^^^ 
let {ut,vt) be a solution of dSBM(^, 

f^)uo,vo CLnd {^ui,Vi) be a solution of 

dSBM(^ 

i'^)ito,iJO' Then the following holds: 

= E"0'^'o [H{uo + vo, UQ - vo, ut + it,ut- Vt)]. 
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Analogously, the self-duality relation holds for the nonspatial model with 
duality function 

H^{u,v,u,v) = exp(— y^l — quu + + Qvv)^ 

mapping (]R>o x ]R>o)^ to C. 

4. Weak longtime convergence. In this section we discuss weak longtime 
convergence of symbiotic branching models and prove Proposition 2.1. We 
proceed in two steps: first, we prove convergence in law to some limit law 
following the proof of [7] for = 0. Second, to characterize the limit law for 
the spatial models, we reduce the problem to the nonspatial model. 

Proposition 4.1. Let q^{-\,X),k>^ and{ut,Vi) a solution of either 
cSBM(e , k)uq,dq or dSBM(g, K)uo,t;o with initial conditions uq — u,fo — v. 
Then, as t^ oo, the law of {ut,vt) converges weakly on M^^^ to some limit 

(Uqo , f oo ) • 

Proof. The proof is only given for the discrete spatial case and the 
continuous case is completely analogous. Let us first recall the strategy of 
[7] for Q = which can also be applied with the generalized self-duality re- 
quired here. Convergence of {ut,vt) in M^^^ follows from convergence of 
{uf + vt,ut — vt) in E. Using Lemma 2.3(c) of [7], it suffices to show conver- 
gence of ¥,^'^[H{ut + Vt,Ut — Vt,(j),ip)] for all € Ef. Furthermore, the 
limit {uoo,Voo) is uniquely determined by E,^'^[H{uoo +Voo,Uoo — Voo,(l),ip)] 
(see Lemma 2.3(b) of [7]). Hence, it suffices to show convergence of 

(4.1) E"'^[i?(ut + Vt, ut - vt,<p, ip)] = E"'^[e-v^<"*+^*''^>+*V^<"*-^*''^>], 

for all {(f),ip) £ Ef. Note that the technical condition of Lemma 2.3(c) of [7] 
is fuUfilled since due to Proposition 3.1 

E^'^[{ut + vt,(l)-x)] = {u + v){l, Pt^.x) < C < oo. 

To ensure convergence of (4.1) we employ the generalized Mytnik self-duality 
of Lemma 3.5 with uq := ^^-^,Vo := 

jgu.v ^^-^l-g{ut+vt,4>)+i^/l+g{ut-vt,tl!)j 

— jgwojUo ^^~y'l~g{ut+vt,uo+i>o)+iV'^+Q{'^t-vt,uo-vo)j 

(4.2) _ _ 

— jgitOj-Co ^g-Vl-e("+f)(l,Mt+t't>+iv'l+£'("-f)(l:"t-i't> j 

By assumption, no,'Uo have compact support and hence by Proposition 3.2 
the total-mass processes {l,ut) and (l,^*) are nonnegative martingales. 
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By the martingale convergence theorem {l,ut) and (1,^4) converge almost 
surely to finite limits denoted by (1,-Uoo), {l,Voo). Finally, the dominated 
convergence theorem implies convergence of the right-hand side of (4.2) to 

(4 3) ]£"0,«'0 Jg-v'l-£'("+'-')(l:'"oo+'i'cx)>+Vl + £'(«-'!')(l>"cx)-Soo>j 

Combining the above, we have proved convergence of 
jgu, V ^g- vi-e(«t +f t ,<l>) +i\/i+g{ut-vt ,-!/)> j 

which ensures weak convergence of {ut,vt) in M^^,^ to some limit which is 
uniquely determined by (4.3). □ 

Again, as in Remark 2.2, the previous proposition can be proved for non- 
deterministic initial conditions as in [4]. 

The rest of this section is devoted to identifying the limit (uoo,t'oo) in 
the recurrent case. Before completing the proof of Theorem 2.1 we discuss 
a version of Knight's extension of the Dubins-Schwarz theorem (see [14], 
3.4.16) for nonorthogonal continuous local martingales. 

Lemma 4.2. Let (Nt) and {Mt) he continuous local martingales with 
Nq = A-Iq = almost surely. Assume further that, for t>0, 

[M.,M.]t = [N.,N.]t and [M. , N.]t = q[M. , M.]t a.s., 

where g£ [—1,1]. If [M. , M.]oo = 00 a.s., then 

is a pair of Brownian motions with covariances [B^,B'^]t = gt, where 
(4.4) T{t)=mi{s:[M.,M.],>t]. 

Proof. It follows from the Dubins-Schwarz theorem that B^,B^ are 
each Brownian motions. Further, by the definition of T{t) we obtain the 
claim 

[B^ , B% = [M. , N.]T(t) = q[M. , M.]Tit) = Qt. □ 

Remark 4.3. If T* := [M,M.]oo < 00 the situation becomes slightly 
more delicate but one can use a local version of Lemma 4.2. Indeed, define, 
for t > 0, 

^ ^ *■ \Mt*, fort>r*, 

where the time-change T is given in (4.4) and define B^ analogously for N 
(recall that [M.,M.]t = [N.,N.]t). Then the processes B^,B^ are Brownian 
motions stopped at time T* . The covariance is again given by 

[B\B%;,T*=g{tAT*), t>0. 
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For the rest of this section let B^,B^ be standard Brownian motions with 
covariance 

(4.6) [B\B% = gt 

started in u,v, denote their expectations by i?"'^', and let 

T = M{t:B}B^ = 0}. 

The above discussion can now be used to understand the longtime be- 
havior of symbiotic branching processes. We start by giving a proof for the 
nonspatial symbiotic branching model and then modify the proof to capture 
the corresponding result for the spatial models. 

Proposition 4.4. Let {ut,vt) be a solution of SBM{q, k)u,v Then, as 
t — )• oo, {ut,vt) converges almost surely to some (uqc'Voo)- Furthermore, 
C^''{uoo,Voo)=iy'''{Bl,Bl) with Bl,Bl from Proposition 2.1. 

Proof. Solutions of the nonspatial symbiotic branching model are non- 
negative martingales and hence converge almost surely. This implies the first 
part of the claim and it only remains to characterize the limit. Obviously, the 
L^-martingales {ut), [vt] satisfy the cross- variation structure assumptions of 
Lemma 4.2 and, thus, {ut,vt) = {B^_i^^-^, B^_i^^^). To obtain the result, we 

need to check that r~^(oo) = r. By definition of SBM, the time-change is 
given by 

(4.7) T-\t) = [u.,u.]t = 

To see that T^^(oo) = r < oo, first note that T~^(i) < r for all t>0. This 
is true since ut = B^_^^_^yVt = B'^_i^^^ and solutions of SBM are nonnega- 

tive. To argue that T~^{t) increases to r, more care is needed. Since the 
martingales converge almost surely, T~^(t) converges to some value a <t. 
Suppose a <T, then {ut,vt) converges to some {x,y) with x,y > 0. This 
yields a contradiction since T~^{t) = k f^Usygds would increase to infinity. 
Hence, almost surely, 

{ut,vt) = (5^-i(j),-b|_i(^)) *^ (-b^-i(oo),-b|-i(oo)) = {bI.bI). □ 

In particular, the proof of Proposition 4.4 provides an important relation 
for {BI,B^). As remarked below Lemma 3.5, the self-duality also works in 
the nonspatial model: 

]g"0 '^0 ^g- ^/l-Q{ut +vt ) {uo+vo)+iVi+d{ut-vt ) (So --Co ) j 

_ wuo,% \„-V^-(>i'^o+'Uo){ut+vt)+i^l+g(uo-vo){ut-vt)] 



^KUgUs dB^, / y/KUgVs dB 



K. / UgVs ds. 
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where both {ut,vt) and {ut,vt) are solutions of SBM(f3, k) with different 
initial conditions. As shown in the proof of Proposition 4.4, {ut,vt) [resp., 
{ut,vt)] converges almost surely to {B^^B'^) with initial condition {uo,vq) 
[resp., {uo,vo)]. Using dominated convergence, this shows the following du- 
ality relation for (B^^B^) when started in initial conditions {u,v), {u,v): 

E^'^[H^{Bl + B^, Bl - B^,u + v,u - v)] 

(4.8) 

= E^''' [H\Bl + Bl,Bl-Blu + v,u- v)]. 

Proof of Proposition 2.1. Again, the proof is only presented in the 
discrete spatial setting since the continuous case is analogous. We retain the 
notation of the proof of Proposition 4.1 where we showed that, as t tends to 
infinity, 

jgu,v j-g-v'l-£'("t +''*></')+« v'l+£'('f j 

Let us specify the limit law as for the nonspatial symbiotic branching pro- 
cess. As seen in Proposition 3.2 the total-mass processes ut := {ut,V) and 
vt '■= {vt, 1) are nonnegative continuous L^-martingales with cross- variations 
= g[u.,u.]t = Q[v-,v.]t,t > 0. Thus, by Lemma 4.2, reasoning as in 
(4.7), {ut,vt) = iB^^i^_i,y B^_i^^^), where B^,B'^ are Brownian motions started 

in Uo = {^-^, 1), vo = (^^) 1) with covariance [B^,B^]t = gt and T~^{t) = 
K ^^{ustVs) ds. Again, we need to show that T~^(oo) = r. This is much more 
subtle than in the nonspatial case since the quadratic variation might level 
off even if both total- mass processes ut, vt are strictly positive. Li [7] it was 
shown that for Q = 0, almost surely, this does not happen in the recurrent 
case [cf. the proof of their Theorem 1.2(b)]. Their proof can be used directly 
for q£ ( — 1,1). Hence, almost surely, 

(4.9) {{nt,l),{vul)y^{Bl,B'^). 

Combining the above discussion with (4.3), we are able to determine the 
limit. First, we derived 

^((<^+'/')/2,l>,({?i-^)/2,l) jg- v^(«+i))(Bi +B?)+iv^(«-")(Br--B?)] ^ 

To use Lemma 2.3(c) of [7] we manipulate the right-hand side using (4.8): 

_g((<i!.+t/')/2,l>,((</,-V.)/2,l)jg-v^(«+t))(Bi+B?)+iv^{«-«){B^-S?)] 

^ ^((<^+V.)/2,i>,((</,-V>)/2,i> [^0 ^Bl + BlBl-Blu + v,u- v)] 
= E^'^[H\Bl + Bl Bl - Bl 1), ii^, 1))] 
= E-'^[H{Bl + BlBl-Bl<p,^P)], 
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where, as in Proposition 2.1, (resp., B^) denotes the constant function 
taking only the (random) value B^ (resp., B'^). In total we have 

K^'-'iHiut + vuUt - vt,<P,n E^^-[H{Bl + BlBl - ^2,<^,^)], 

which implies weak convergence in M"^^^ of {ut,vt) to {Bl,B^) by Lemma 
2.3(c) of [7]. □ 

5. Moments. In this section we prove Theorems 2.5 and 2.7. Before giv- 
ing the proofs we prove an equivalence for moments of correlated Brownian 
motions. 

5.1. Moments of the exit-point and exit-time distribution of correlated 
Brownian motions in a quadrant. Let g G (—1,1), u,v >0 and B^,B^ be 
Brownian motions started in u,v with 

(5.1) {B\B^)t = Qt. 

The starting points of Brownian motions will be indicated by superscripts 
in probabilities and expectations. Further, let 

(5.2) =mi{t>0:BlB^ = 0}. 



Theorem 5.1. Let p > and u,v > 0. Under the above assumptions, 
the following conditions are equivalent: 

(i) 

TT 

7r/2 + arctan(£i/(^l - £>2)) ' 

(ii) 
(iii) 

E-'^[\{Bl,,Bl,)\P]<^. 

Proof. We start with the proof of the equivalence of (i) and (ii). Define 
a cone in the plane with angle 9 G (0, 2tt) by 

C{ip) = {re''^ : r > 0, < (/) < 

and denote its boundary by dC{ip). Note that with this definition, the posi- 
tive real line is always contained in C{ip). Further, we define, for g £ (—1,1), 
a sector in by 

S{q) = |(x,y) G .X > 0,y > --^J==x^ 
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and denote by dS{g) its boundary. Note that this time, the positive imagi- 
nary axis is always in S{q) and that the angle of the sector at the origin is 
given by 



TT 

— h arctan 
2 



(^/T^)' 



To transform the correlated Brownian motions B^,B^ to planar Brown- 
ian motion we use the simple fact that := B^,W^ := ) defines 

a pair of independent Brownian motions started in u, ( ) satisfying 



{B^,B'^) = {W^,qW^ + yJl-Q^W'^). By the definition of S{q), the planar 

Brownian motion started in iu, ( '"zB^ )) hits dSig) if and only 

yi-r 

if the correlated Brownian motions B^,B^ started in u,v hit dC{^). Hence, 
for as in (5.2), we have 

(5.3) = := inf{t > : {W^,W^) G dS{g)}. 

Since planar Brownian motion is rotation invariant, S{q) may be rotated to 
agree with the cone 0(0), without changing the exit time. Obviously, with 
the corresponding rotated initial conditions, the law of the first exit time 
Tc(^0^ from the cone C{9) agrees with the law of r^. For planar Brownian 
motion in a cone C{0) it is well known (see [23], Theorem 2) that 



(5.4) E-'y[{Tc^g)f'] <oc ^ p< 



vr 



independently of x,y. (5.3) and (5.4) now imply the equivalence of (i) and 
(ii) and independence oi u,v. 

The proof of the equivalence of (i) and (iii) is via conformal transformation 
of the cone C{9) to the upper half-plane. Indeed, we are going to calculate 
the densities of the exit-point distributions 

(5.5) P^^^{Bl^=0,B^^^>y), P"'"(i?^s < x, i?^^ = 0). 

We proceed in three steps: after reducing to independent Brownian motions 
in S{g) as for the exit time, we rotate S{g) to C{9) and, finally, stretch the 
cone to end up with the upper half-plane. 

Recall that the first exit of {B^,B'^) happens at position (0,y) G <9C(^) 
if and only if the first exit of {W^,W'^) takes place at (0, —t^ — ) G dS{g). 

Hence, (5.5) transforms to 

p-^-{Bl, = 0,Bls>y) 

(5.6) 
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In a similar fashion one obtains 

— p",{^'-£'«)/\/i-e^ 



(5.7) 



< 



We represent the transformed initial conditions (^1,2:2) = {u, 
in polar coordinates, that is, 



V — QU 



Z2 



/ „ iv- QuY ( 

' H K — cos arctan 

1- V 



( 



QU 



{v — guy 



sin ^arctan ^ 



- QU 



For the rotation we add the angle arctan ( 



to get the new initial 



condition. Finally, to map the cone C{9) conformally to the upper half-plane 
H, we apply the map z 1— t- z'^^^ which maps C{6) onto H. Using conformal 
invariance of Brownian motion (e.g.. Lemma 7.19 of [15]), the problem is 
reduced to the computation of the exit distribution of planar (time-changed) 
Brownian motion from the upper half-plane. Indeed, due to the random time 
change the (almost surely finite) exit time changes but not the distribution of 
the exit points, which is Cauchy (see Theorem 2.37 of [15]). Thus, to obtain 
the distribution of the exit points explicitly it, only remains to specify the 
transformed initial condition zi,Z2, which is given by 



n2 + 



Z2 



X COS ( — I arctan 



V — QU 



Q^U 



+ arctan 



X sm 



tt( ( 
— arctan 



V — QU 



+ arctan 



Now, let W^,W'^ be two independent Brownian motions with Wq = zi,Wq 
Z2 and 



r^:=inf{t>0:Wi^ = 0}. 



ON THE SYMBIOTIC BRANCHING MODEL 25 
Then, by (5.6), (5.7), 

P'^'^'iBl^ = 0,i?2^ >y) = P^'(-~(^^yV^(w^^: = 0, W^^ > ^ y ] 

7r/0N 



pzi,~z2l o<T^\,, < ( x( 1 + 



Q 



2 \ l/2\ TT 



l-£»2 

n/e 



Explicit manipulations of the Cauchy distribution yield 
(5.8) 

1 vr/eW^-i 

ar, 



y ^z2A/r^"^ l + (((^/Vl^)'^/' + ^i)/^2)2 

^"'"(5^5 <X,B2^ =0) 

(5.9) 

1 7r/6lW^-i 



■ dr. 



Finally, noting that Jq°° '^j^i^"^^^ da; < oo if and only if p < a, we deduce from 



that 

(5.8) and (5.9) that 



E"'''[|(5is,5^s)r]<oo if and only if p<^. □ 

5.2. Proof of Theorem 2.5. The proof relies on a combination of the 
self-duality based technique of the proof of Proposition 2.3 and the close 
relation between the moments of the exit-time and exit-point distribution 
of correlated Brownian motions obtained in Theorem 5.1. 



Proof of Theorem 2.5. We proceed in several steps. First, the result 
for the nonspatial model is proved and thereafter the results for the discrete- 
space and the continuous-space models. Finally, we present the argument in 
the transient case. In the following we use the definition of B^,B^ and r 
from Proposition 2.1. 

Step 1. Suppose {ut,vt) is a solution of SBM(^, k)i^i. 
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"=^": We first assume g < q{p), in which case Theorem 5.1 imphes 

< oo. As argued in the proof of Proposition 4.4, ut is a non- 
negative martingale and due to the same arguments satisfies E^'^[[n.]^''^] < 

< OO for all i > and k > 0. Considering ut = ut — UQ = ut — 1, we 
apply the Burkholder-Davis-Gundy inequality to get 

Ei'i[<] = Ei'i[(uj + l)P] 

= Ei'i[l|^^<i|(^2i + If] + Ei'i[l{^,>i|(nt + If] 



<Cy + C„E 



1,1 



sup vFg 



i-0<s<t 

<Cp + C'^¥}'\[u.]l'^]<oo 

independently of t and k. 

"<;=." Conversely, for g > q{p), Theorem 5.1 implies that E^'^[{BlY] = oo. 
Using Fatou's lemma and almost sure convergence of ut to B^, the proof for 
the nonspatial case is finished with 

liminf E^'^uf] > E^'^u*^] = E^'^[{Brf] = oo. 

Again, this lower bound is independent of k. 

Step 2. The proof for dSBM(^>, is started by reducing the moments 
for homogeneous initial conditions to finite initial conditions. Indeed, em- 
ploying Lemma 3.5 with = ^ = where 1^ denotes the indicator func- 
tion of site A; G Z*^, gives 

— ^g-v'l-£'(l + l:"t+^'t>] 

where we used the argument of Remark 2.6. Note that, due to our choice 
of initial conditions, the complex part of the self-duality vanishes. Since the 
above is a Laplace transform identity, we have 

C^'Hut{k)+vt{k))=/:^'"^>'{{l,ut) + {l,vt)) 

and hence 

(5.10) Ei'M(^zt(fe) + i;t(A:)n=Ei'='i'=[((l,nt) + (l,^i)n. 

We are now prepared to finish the proof of the theorem for the discrete case. 

Suppose g < g{p). Let Mt = {l,ut) + (l,^*), which due to Lemma 
3.2 is a square-integrable martingale with quadratic variation 



[M.]t = [{l,u.)]t + [{l,v.)]t + 2[{l,u.),{l,v.)]t = (2 + 2g)[{l,u.) 
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To apply the Burkholder-Davis-Gundy inequality, we switch again from M 
to Mt = Mt — Mq, which is a martingale null at zero. Hence, 



E 



[Mf 



[{Mt + Mof ] <Cp + CX'' ' [Mf] 



Then we get from (5.10) and the Burkholder-Davis-Gundy inequality 



<Cp + CpE^*'^^- 



sup Mf 

0<s<t 



for some constants Cp, Cp independent of t and k. As in the proof of Theorem 
2.1, the random time change which makes the pair of total masses a pair of 
correlated Brownian motions is bounded by r, that is, [{l,u.)]t < r for all 
t > 0. This yields by Theorem 5.1 

E^'^[utik)P] < E^^^iiutik) + vt{k))P] <Cp + C'p{2 + 2Qy/^E^^\T'p'^] < oo. 

"<;=." Suppose Q > q{p). As in the proof of Theorem 2.1 we use the almost 
sure convergence of ((1, ■Uf), (1, vt)) to {Bl,B'^). Combining this with Fatou's 
lemma gives 

liminfE^^-'i''-[((l, nt) + (1,^*))^] > liminf E^'-^'^ [(1, 

t—>-oo t-^oo 



liminf (1, ttf)^ 

t~>-co 



The right-hand side is infinite due to Theorem 5.1 and hence E^'='^'= [((1, ut) + 
{l,vt)y] diverges. Equation (5.10) now shows that E^'^[{ut{k) + vt{k))P] also 
grows without bound. Since symbiotic branching processes are nonnegative, 
this is also true for E^'-^[ut{ky] as can be seen as follows: 

Ei'i[(^.i(A;) + z;t(A;))T<Ei'M(2nt(A;))n|„,(,)>,,(,)|] 

+ E^'^[{2vt{k))n{^,(^k)<vm}] 

< 2PE'^^^[ut{k)P] + 2PE^^^ [vt{k)P] 

= 2P+^E'^'^[utik)P], 

where we used Lemma 3.3 to see that E^''^[ut{ky] =E^''^[vt{k)'P]. 

Step 3. The proof for cSBM(^, k)i^i is slightly more involved since we 
cannot use the indicator Ix to get ut{x) = {ut,lx), where now {f,g) = 
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J^f{x)g{x)dx. Instead we use a standard smoothing procedure. For fixed 

j; e M let 



where we skip the dependence on x. The main part is to show that 

^^{Ut{x) +Vt{x)) - {{ut,Pe) + {vt,Pe))\\LP 



(5.11) 



which imphes 



< \\ut{x) - {ut,Pe)\\LP + \\vt{x) - {vt,Pe)\\LP 0, 



(5.12) 



lim ||(«t,Pe) + {vt,Pe)\\LP = \\ut{x) +Vt{x)\\LP■ 
s—^0 



Due to symmetry we only consider ||nt(x) — (nt,Pe)||LP. To prove (5.11) 
we first observe that, due to the Green function representation provided in 
Corollary 19 of [8], 

\\ut{x) - {ut,Pe)\\LP 

PtUo{x) - {Pt+e,UQ) + [ [ Pt-s{x - b)M{ds,db) 

Jo JR 



Pt^sPeix - b)M{ds,db) 



LP 



where M(ds, db) is a zero-mean martingale measure with quadratic variation 




f{s,b)M{ds,db) 



f'^{s, b)us{b)vs{b) ds db 



for test functions / such that the integral is well defined (see Lemma 18 of 
[8] for details). 

For homogeneous initial conditions, the first difference vanishes and it 
suffices to concentrate on the difference of the stochastic integrals. By the 
Burkholder-Davis-Gundy inequality the difference of the integrals can be 
estimated as 



E 



Pt-s{x - b)M{ds,db) 



Pt-sPe{x-b)M{ds,db) 



{pt-s{x - b) -pe+t-six - b)) Us{b)vs{b)dsdb 



p/2' 



Now expanding {pt-s{x - b) -pe+t-s{x - b))'^Us{b)vs{b) as 
{pt-s{x-b)-pe+t-s{x-b)f^-^'^'^ 

X {pt-six - b) -pe+t-s{x - b)f^'^Us{b)Vs{b), 



ON THE SYMBIOTIC BRANCHING MODEL 29 

we get the upper bound (taking the expectation under the integral is vahd 
since the integrands are nonnegative) 

p-i 



{pt-s{x - b) -pe+t-s{x - b)fdsdb 



{pt^six - b)-pe+t-s{x - b)fn{M^>s{W] dsdb 



where we have used that, for f,g£LP, 

(^2(p-i)/p)(^2/p^)^^y < (1 fdxj 'I fgPdx 

by Holder's inequality. As in [8], page 153, the second term can now be 
bounded from above by a constant depending only on p and t. The first 
factor can be estimated by e(P~^)/2 due to [22], Lemma 6.2. Hence, for fixed 
p > l,x gM and t > 0, (5.11) holds and thus we obtain (5.12). The rest of 
the proof is similar to the discrete case but slightly more technical. Since 
Ps{x — ■) is rapidly decreasing, we have 

Thus, we get 

C^'\{nt + vt,Pe)) = C^^'^^{{l,nt + vt)) 

and in particular 

E^'^iiiut + vt,Pe)r]=w^'P^[{{i,ut) + {i,vt))n- 

We may now finish the proof in a similar way to the discrete case. 

"=^>." Due to (5.12) we are done if we can bound K^''^[{ut + in- 
dependently of e > and t > 0. This can be done as before: {l,ut) and 
{l,vt) are random time-changed correlated Brownian motions with initial 
conditions {l,Pe) = 1 for all e > 0. Using, as before, the auxiliary martingale 

Mt = {l,ut) + {l,vt) - (l,uo> - {l,vo), 

we obtain (as in the discrete case) with the help of the Burkholder-Davis- 
Gundy inequality 

E^^^utix) + vtix))P] = hm E^^^iiut + vt,Ps)''] 

£->0 

<Cp + C7plimEP-P-[Mf] 
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<Cp + C7;(2 + 2^f/2^i.i[^p/2]. 

The positive constants Cp, C'p are independent of e and t, whereas M and 
the random time change [M.]t do depend on e. However, the bound [M.]t < t 
holds for all e > and t>0 since Bq = Bq = {l,pir) = 1. For g < q{p) the 
right-hand side is finite by Theorem 5.1 and independent of t > 0. Since 
E-'-'-'-[nt(2;)P] <E,^'^[(ut{x) + vt{x))P], the first direction is shown. 

"<;=." First note that by translation invariance of initial condition, spatial 
motion and white noise 

E'^^[{ut{x) + vt{x)r] = Ei'i[(nt(y) + vt{y)r] 

for fixed time t > and arbitrary spatial positions x,y G M implying that 

rx+l/2 

E^'\ut{x)+vt{x)y]= / E^'^[{ut{y)+vt{y)y]dy. 

Jx~l/2 

Using Fubini's theorem and Jensen's inequality we obtain for p> 1 the lower 
bound 

fX+l/2 \ p- 



Ei'i[(nt(x) + t;i(x)n>Ei'i 



{ut{y) + vt{y)) dy 



We now choose an arbitrary nonnegative (nontrivial) smooth function / with 
support contained in {x — l/2,x + 1/2) that is bounded by 1 and integrates 
to some c G (0, 1), say. A lower bound is now given by 

E^'^[{ut{x) + vt{x)y] > E^'^iint + vt, fY] 

= Ef'f[{ut + it,ir] 
>E/'/[(nt,in, 

where we utilized for the equality the self-duality of Proposition 5 of [8] . 

Finally, as in the discrete case, Fatou's lemma and the martingale con- 
vergence theorem imply 

\immfE^'^[{ut{x)+Vt{x)f] > ] = oo 

by Theorem 5.1 and due to nonnegativity of solutions as well 

liminfE^'^[ut(x)P] = oo 

t—>-co 

proving the claim. 

Step 4. The first direction of the above proof for dSBM(£), also works 

for the transient case since Eifc'ifc[[Af.]^^] <£;i'i[rP/2] is independent of 
recurrence/transience. □ 



ON THE SYMBIOTIC BRANCHING MODEL 



31 



5.3. Proof of Theorem 2.7. We now study the "criticality" of the crit- 
ical curve in more detail. As a preliminary result (mixed) moments of the 
nonspatial model are analyzed. The idea is to combine three different tech- 
niques: the martingale argument which led to Theorem 2.5 for E^'^ftiJ'], a 
perturbation argument based on the moment duality which allows us to 
deduce exponential increase/decrease of E^'^[it"~^ft], and finally moment 
equations which yield exponential increase/decrease for all mixed moments 

Proposition 5.2. The following hold for nonspatial symbiotic branch- 
ing processes: 

(1) For all K>0 and n G N.- 

• E^'^[u"] grows to a finite constant if g < g{n), 

• E^'^[u"] grows subexponentially fast to infinity if g = g{n), 

• E^'^[uj'] grows exponentially fast if g> g{n). 

(2) For all K>0, n S N and m = I, . . . ,n - 1: 

• E^'-'^[ii"~™f™] decreases exponentially fast if g < g{n) , 

• K^'^[u^~"^v^^] neither grows exponentially fast nor decreases exponen- 
tially fast if g = g{n), 

• E^'-'^[ii"~'"f™] grows exponentially fast if g > g{n) . 

Proof. Step 1. Martingale arguments based on the connection of mo- 
ments of exit times and exit points of correlated Brownian motions were 
carried out in the proof of Theorem 2.5. This led to the first part of (1). 
Applying Holder's inequality with p = , q= we get the bound 

(5.13) E^'i[<-"*wn < E^'iK]("-™)/"E^'iK]™> = E^'^uJ^] 

by symmetry. This implies that for g < g[n) all mixed moments stay bounded 
as well. 

Step 2. We apply the moment duality for the nonspatial model as ex- 
plained in Remark 3.4. Combining the duality with the martingale argument 
of the first step we can understand the case g < g{n) for mixed moments 
in a simple way. Note that for mixed moments the dual process starts with 
n — m particles of one color and m particles of the other color at time 0. 
Note that for mixed moments Lf > t, since there is always at least one pair 
of different color. Now suppose g < g{n), then for < e < g{n) — g we get 

Ei'i[u^-'"t;™] =E[e^(-^^+^'^t'')] =E[e^(-^^+(^'+=)^^)e^^=-^^] 

< E[e'^'^^^~^^^~^'^^^t''^]e~^'^K 

Since the first factor of the right-hand side is just the moment E,^'^[u^~'^v^] 
for g + £ strictly smaller than g{n), this is bounded for all t and k. Hence, 
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for Q < Q{n) all mixed moments decrease exponentially fast proving the first 
part of (2). Note that since is a submartingale, the moment E^'"'^[u"] is 
nondecr easing. 

For Q = Q{n) we first consider the pure moments. Again, for the critical 
case, Theorem 2.5 implies 

j£[g«(Lr+(e(n)-£)Lf )] < < ^ 

for all e > and t > 0. With the crude estimate hf < (^^'^t we get 

Since e is arbitrary this implies subexponential growth to infinity of E^'^ [u^ (/c)"] 
at the critical point. Hence, the second part of (1) is proven and combined 
with (5.13) so is the upper bound of the second part of (2). 

Step 3. A direct application of Ito's lemma and Fubini's theorem yields 

E^'^K] = 1 + K 2 ) J^^^'^K~^ys] ds. 

Since we already know from the martingale arguments that E^'"'^[u"] increases 
to infinity in the critical case, the mixed moment E[ti"~^z;t] cannot decrease 
exponentially fast proving the lower bound of part two of (2). Furthermore, 
with the same arguments as above, for g > Q{n), this leads to 

Since the first factor of the right-hand side equals E[u"~^t;f] at the critical 
point, it does not decrease exponentially fast. Hence, the product increases 
exponentially fast. In particular, due to (5.13), this also implies the third 
part of (1). Now it only remains to prove exponential increase for the other 
mixed moments. Again, using Ito's lemma and Fubini's theorem yields the 
following moment equations for the mixed moments: 

E^'^[u'^~\^] = l + K fE^'^[u'^~\s]ds + g{n-2)K f E^^\vJl-^vl]ds 
Jo Jo 

and similarly for all other mixed moments. Since we already know that 
E^'^[ti"~^t;f] grows exponentially fast in i, this implies exponential growth of 
E^'^[it"~^fj]. Iterating this argument gives exponential growth of all mixed 
moments for g > g{n). This shows the third part of (2) and the proof is 
finished. □ 



Now it only remains to prove Theorem 2.7, where some ideas for the 
nonspatial case are recycled. 



ON THE SYMBIOTIC BRANCHING MODEL 



33 



Proof of Theorem 2.7. First, due to Lemma 3.3, for homogeneous 
initial conditions, the moments of ut{k) and vt{k) are equal for all t > 0. For 
the existence of the Lyapunov exponents we use a standard subadditivity 
argument. Hence, it suffices to show 

Using Lemma 3.3, we reduce the problem to E[e''(^i +£'i'f )]^ where the dual 
process (rit) starts with n particles of the same color all placed at site k. By 
the tower property and the strong Markov property, we obtain 

We are done if we can show that 

(5.14) E^'je'^C-^-r+eif )j < jEwojgKCir+e-f'f )] 

for any given initial configuration n' of the dual process consisting of n parti- 
cles. The general initial conditions of the dual process consist of particles 
of one color and particles of the other color (n^ = n) distributed arbi- 
trarily in space at positions ki, . . . ,kn- Using the duality relation of Lemma 
3.3, we obtain 

jgn'jg.(Lr+,Lr)] ^ E^'^[Us{ki) ■ ■ ■ Us{k^i)Vs{k^i+i) ■ ■ ■ (fe„l+„2 )] 

where, in the penultimate step, we have used the generalized Holder inequal- 
ity. 

Having established existence of the Lyapunov exponents, we now turn to 
the more interesting question of positivity. The boundedness for g < Q{n) in 
Theorem 2.5 immediately implies that in this case ^{g, k) = 0. Now suppose 
Q= Q{n), that is, {Q,n) lies on critical curve. We use the perturbation argu- 
ment which we already used for the nonspatial case combined with Lemma 
3.3 and Theorem 2.5 to prove that in this case moments only grow subexpo- 
nentially fast. This implies that the Lyapunov exponents are zero. Again we 
switch from E^'^[ut{k)''] to E[e'^^^^+^^t -*], where the dual process is started 
with all particles at the same site and the same color. Since moments below 
the critical curve are bounded, we can proceed as for the nonspatial model. 
For any e > 0, we get 

oo > C{e) > Ele'^^^^+^^^^e-"'^^] > Eie^'^^^+^^^^e-^'^"^)* 

>E'^'^[ut{kr]e-''<2)\ 

where we estimated the collision time of particles of different colors by the 
collision time of all particles which is bounded from above by (2)^- Since e 
on the right-hand side is arbitrary, 7(^, k) cannot be positive. 
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Finally, we assume g > Q{n). The idea is to reduce the problem to the 
nonspatial case which we already discussed in Proposition 5.2. Actually, 
we prove more than stated in the theorem since we also show that mixed 
moments E-*-'-*- [ut{k)'^~'^vt{k)"^] grow exponentially fast. For m = 1, . . . , n — 1 
the perturbation argument leads to 

¥}''^[ut{kT-''^vt{k)"'] = E[e''(^^+^^^)] = ^[e''^Lf+e{n)Lf)^K(e-e{n))Lf^^^ 

The idea is to obtain a lower bound by conditioning on the event that all 
particles have not changed their spatial positions before time t (but, of 
course, have changed their colors). Under this condition the particle dual is 
precisely the particle dual of the nonspatial model. More precisely, we get 
the lower bound 

^\(,<Lf+g{n)Lf)^K{e-Q{n))Lf.^^ Spatial change of particles before time t] 

= ]E[e'^(^r+f{n)if)g«:(e-e(n))if |no spatial change of particles before time t] 

X P[no spatial change of particles before time t] 

= W.[e^iLf+e{n)Lf)^K{e-Q{n))Lf^^^^ spatial change of particles] e""*, 

where the final equality is valid since the event {no spatial change of particles 
before time t} has probability e~"*. This is true since the event is precisely 
the event that n independent exponential clocks with parameter 1 did not 
ring before time t. For 1 < m < n — 1 there is always at least one pair of 
particles of different colors and, hence, we get the lower bound 

]gjgK{Lr+e(n)Lf)|^Q spatial change of particles until time t]e'^(^^-''("))*e-"*, 

which equals 

E^'i [u'l-"'v'^]e''^s-e{n))t^-nt 

for a nonspatial symbiotic branching process with critical correlation g = 
Q{n). Choosing k such that k{q — Q{n)) > n the result now follows from 
Proposition 5.2. □ 

As mentioned in the course of the proof, we actually proved that for 
Q > g{n) and m = 0, . . . , n 

E^'^[ut{kr-"'ik)vr{k)] 

grows exponentially in t. As for the nonspatial model one could ask whether, 
and if so how fast, mixed moments decrease for g < g{n). For the second 
moments it was shown in [1] that for g < g{2) = 

71; 

Ll, ci>3, 



E^^^[ut{k)vtik)] 
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where fa denotes weak asymptotic equivalence as f — )• oo. It would be inter- 
esting to see whether or not different rates of decrease appear for moments. 

A detailed quantitative study of the Lyapunov exponents as functions of 
Q and K has so far only been carried out for second moments (see [1]). In 
contrast to the parabolic Anderson model, where higher Lyapunov exponents 
are well studied (see [9]), we do not have much insight. Only a first upper 
bound for the Lyapunov exponents in n and the distance to the critical curve 
can be obtained from the perturbation argument of the previous proof. 

Proposition 5.3. IfQ>g{n), then -fn{Q,K) < K{^){g - g{n)). 

Proof. By Lemma 3.3 and Theorem 2.5 for g > g{n), there are con- 
stants C(e) such that 



yielding the result. □ 

6. Speed of propagation of the interface. In this section we show how to 
use the moment bounds of Theorem 2.5 to obtain an improved upper bound 
on the speed of propagation of the interface as defined in Definition 2.10. 
We will only sketch the crucial parts in the proof of Theorem 6 of [8] that 
need modification. Note that the method used here is based on Mueller's 
"dyadic grid technique" introduced in [17]. 

Proof of Theorem 2.11. To prove that the interface wih eventually 
be contained in 



(for suitable C > 0), by symmetry, it suffices to show that the right endpoint 
of the interface 




Hence, for all e > 



Ei'i[nt(A;)"] < C(e)e''(^-^(")+"K2)* 



[-Cv/riog(r), cVriog(r)] 



R{ut) := sup{3; e M|nt(x) > 0} 



up to time T can eventually be bounded by C\jT log(T). To this end we 
define 
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and show that, for suitably chosen C, P^it"'^K+ (hmsup^gp} A^) =0. By the 
Borel-Cantelh lemma, this follows from 

oo 

(6.1) ^pV-.1m+(A„) <oo. 

n=0 

In the following we modify the arguments of [8] to obtain an upper bound 
for F^v-^\+(^An) which is sumamble over n. 

Lemma 6.1. For any integer n there is a finite constant Cn such that for 
Q < g{4n - 1) 

E^K- =1*+ [{ut{x)vt{x))"] < CnVPdR-{x), xeR,t>0. 

Proof. First recall from (87) of [8] that E^k- [ut{x)] = Ptlj^- (x). We 
now use Holder's inequality and Theorem 2.5 to reduce the mixed moment 
to the first moment: 

EV-.V+[(ut(3;)t;i(x))"] 

= EV-ia+ [ut(2;)^/2^j(x)"^i/2^t(x)"] 

< (e1k-'1k+ [nt(x)])i/2(Ei'i[nt(x)*"-i])(2"-i)/(8"-2) 

X (Ei'i[t;j(x)^"-i])"/(*"-^). 

This follows from the generalized Holder inequality with exponents 2, {8n — 
2)/(2n — 1) and (4n — l)/n. The first factor yields the heat flow and The- 
orem 2.5 shows that the latter two factors are bounded by constants for 
g<Q{4n-l). □ 

We now strengthen the estimate of Lemma 23 of [8] of the stochastic part 

Nt{b)= f I pt-s{h-a)M{ds,da) 
Jo Jr 

of the convolution representation of solutions of Corollary 20 of [8]. 

Lemma 6.2. For g < g{35) there is a constant C3 such that for e G (0, 1), 
A,T >1, the following estimate holds: 

IPV'1h+ {\Nt{b)\ > e for some t<T andb>A) < Cse"^^— P2t(^)- 

Proof. The proof is along the same lines of [8] replacing only in (116) 
the weaker (exponentially growing) moment bound of [8] by our stronger 
(bounded) moment bound. In the following we sketch the arguments to show 
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where the moments appear. Before performing the "dyadic grid technique," 
increments of Nt need to be estimated. First, by definition 



[\Nt{a) - Nt'{a')\'^'i] 



: E-'-H- '-'-H+ 



{pt^s{b-a) - pt'-s{b - a'))M{ds,db) 



2g- 



which by Burkholder-Davis-Gundy and Holder's inequahty gives the upper 
bound 

[pt^s{b-a) -pt'^s{b- d)f dbds 



[pt-s{b -a)- pt'^sib - o')]'eV.1k+ [(n,(6)r;,(6))''] dbds. 



Ci 



Using Lemma 6.1 and classical heat kernel estimates we can derive (see 
the calculation on pages 153, 154 of [8]) the upper bound 

EiK-'iK+[|iV^(a) -iVi,(a')|^^] 

< C2((|t' - + \a'-a\)A t'/^fiy^tPtl^^ia) + y^t' m^-{a')). 

This upper bound corresponds to (119) of [8] where they have an additional 
exponentially growing factor coming from their moment bound. The dyadic 
grid technique can now be carried out as in [8], choosing q = 9, without car- 
rying along their exponential factor. Hence, we may delete the exponential 
term from their final estimate (110). Note that the necessity of ^> < q{35) 
comes from our choice q = 9 and Lemma 6.1. □ 

The following lemma corresponds to Proposition 24 of [8]. 

Lemma 6.3. If q< i?(35) then, for some constants 0^,0^, the following 
estimate holds for T > 1 and r > C^VT: 



supR{ut)>r)<C5T''pieTir). 



Proof. All we need to do is to argue that Proposition 24 of [8] is valid 
for r > C4VT instead of r > 9^(1 V k)T. We perform the same decomposition 
and note that the estimates of Step 2 of [8] are already given for r > C4VT 
if C4 is large enough. The only trouble occurs in their Step 3. Up to the 
estimate (154), this step works for r > Ci^s/T but here their (weaker) Lemma 
23 produces an exponential in T. More precisely, they need to justify 



T 



22 
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which is only vahd for r > 9^(1 V k)T. As our Lemma 6.2 avoids the expo- 
nential on the left-hand side the estimate holds for r > C4VT with suitably 
chosen C4 and C5. □ 



The significant distinction of the previous lemma to the result of [8] is 
that the inequality is not only valid for r > 9^(1 V k)T but for r > C^VT. 
At this point one might hope to obtain a square-root upper bound for the 
growth of the interface but this fails in the final step in which we validate 
(6.1): 

00 00 

^pV.V+(^„) <^C5n2V6„(C\/nlog(n)) 



n.=0 n=0 
00 



n=0 

00 



22 ^-C^nlog(n)/(32n) 

■v/7r32n 



•-^5 V^„22-C2 



22-C^/32-l/2 



/327r „ 

71=0 



which is finite for C large enough. □ 
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